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The mutual controls of ferroelectricity and magnetism are stepping towards practical applications proposed 
for quite a few promising devices in which multiferroic thin films are involved. Although ferroelectricity 
stemming from specific spiral spin ordering has been reported in highly distorted bulk perovskite 
manganites, the existence of magnetically induced ferroelectricity in the corresponding thin films remains 
an unresolved issue, which unfortunately halts this step. In this work, we report magnetically induced 
electric polarization and its remarkable response to magnetic field (an enhancement of —800% upon a field 
of 2 Tesla at 2 K) in DyMn0 3 thin films grown on Nb-SrTi0 3 substrates. Accompanying with the large 
polarization enhancement, the ferroelectric coercivity corresponding to the magnetic chirality switching 
field is significantly increased. A picture based on coupled multicomponent magnetic structures is proposed 
to understand these features. Moreover, different magnetic anisotropy related to strain-suppressed 
GdFe0 3 -type distortion and Jahn-Teller effect is identified in the films. 



Multiferroics with coexisting and intimately coupled ferroelectric (FE) and magnetic orders, have attracted 
considerable attentions due to their vast application potentials and fundamental interest in the past 
decade 1 " 6 . It has been confirmed that a series of magnetic configurations can break the spatial inversion 
symmetry and thus produce spontaneous polarization (P) 4 . One of the most frequently cited evidences is the 
spiral-spin-ordering (SSO) in orthorhombic manganites £Mn0 3 (R = Tb and Dy), as highlighted by the first 
discovery of multiferroicity and giant magnetoelectric (ME) coupling in bulk TbMn0 3 (TMO) in 2003 7 . For these 
model multiferroics with SSO, successive phase transitions occur due to the proximate free energies of various 
phases in the low temperature (T) region. Upon cooling, the multiferroic £Mn0 3 enters a collinear sinusoidal 
antiferromagnetic (AFM) state at T N (—40 K) from a paramagnetic (PM) one, and then the SSO develops with a 
locked period at T c (—20-30 K) and an improper FE polarization arises simultaneously via the inverse 
Dzyaloshinskii-Moriya (DM) interaction between adjacent spins. The as-generated P can be expressed as P oc 
Cij X (Si X Sj), where denotes the unit vector connecting the two spins S,- and Sj. Upon further decreasing 
temperature to T < 10 K, the moment of rare-earth ions forms a long range ordering at T R and also modulates the 
ferroelectricity 4 ' 8 . 

To date, bulk SSO multiferroic manganites including single crystals have been intensively studied, and the 
associated physics has been progressively understood 9 " 12 . However, research on the thin films is rare, and the 
insights are expected to play a core role in designing advanced multifunctional devices. From the fundamental 
viewpoint, fabrication of thin films will introduce additional degrees of freedom such as strain to mediate the 
multiferroicity besides the bulk-related routes, e.g., chemical doping. It is thus reasonable to expect intriguing 
properties in thin films since the delicate balance among different exchange interactions is sensitive to internal/ 
external perturbations in multiferroic manganites. Indeed, a few recent papers reported emergent ferromagnet- 
ism related to the modified lattice distortion in TMO thin films, while no ferroelectricity associated with the SSO 
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was captured 13 " 18 . This is quite different from the E-type AFM multi- 
ferroic manganite thin films, in which the magnetic ferroelectricity 
was reported in several materials 19 " 21 . 

According to the phase diagram of perovskite manganites, the SSO 
phase is rather sensitive to structure distortion, and can only exist 
within a very narrow region 22 . Furthermore, the SSO in BiFe0 3 was 
revealed to be suppressed by strain 23,24 . It thus seems reasonable to 
expect the suppression of SSO phase by strain in orthomanganite 
thin films. However, while the ground state of bulk YMn0 3 is the E- 
type AFM 20 ' 21,25 , a strain-induced SSO phase accompanied by spon- 
taneous polarization was reported in the corresponding thin films 
recently 26 . In a recent ultrafast optical spectroscopy study, an anom- 
aly around T = 30 K (close to the FE Curie temperature T c ~ 28 K of 
TMO single crystal) was observed in TMO thin films, which hints the 
coexistence of SSO and ferroelectricity but no relevant data were 
reported 27 . These works motivated us to explore the SSO induced 
ferroelectricity and the ME coupling in orthorhombic J RMn0 3 thin 
films. 

Along this line, we choose DyMn0 3 (DMO) thin films as the 
prototype to investigate the possible coexistence of multiple ferroic- 
orders and the ME coupling among them. Besides TMO, DMO is one 
of the typical multiferroic manganites with SSO, and they possess 
quite similar physical properties. Particularly, in both materials, the 
SSO induces ferroelectricity through the inverse DM interaction, and 
magnetic field switches polarization from the c axis to the a axis (P 
flopping) 8 . DMO lies far inside the SSO region in the magnetic phase 
diagram while TMO is somehow in the vicinity of the A-type AFM 
region, allowing more flexibility for spin order modulation in DMO 
than TMO 22 . In addition, DMO single crystal exhibits much larger 
spontaneous polarization (P ~ 2100 |iC/m 2 ) than that of TMO (P ~ 
600 u€/m 2 ) 8 , which facilitates reliable measurements. More impor- 
tantly, in bulk DMO, the Dy moments play a much more active role 
through the Dy-Mn interaction in generating FE polarization than 



the Tb moments in TMO, enabling fascinating physics unavailable in 

TMO 28 " 33 . 

In this work, we perform systematic experiments to investigate 
the structural, magnetic, dielectric, and ferroelectric properties of 
strained DMO thin films grown on Nd-SrTi0 3 (001) (NSTO). 
Magnetically induced ferroelectricity was clearly identified using 
the pyroelectric current method and the direct electric hysteresis 
loop measurement. Interestingly, both the polarization and ferro- 
electric coercivity can be significantly enhanced by magnetic field 
at T < T R , indicating the essential role of Dy moment on ferroelec- 
tricity. In addition, different magnetic anisotropy was observed. 
Notably, these multiferroic behaviors in DMO thin films are quanti- 
tatively different from but qualitatively similar to those for bulk 
DMO crystals, indicating that the strain may play a non-negligible 
role in modulating multiferroicity in such manganite thin films. 

Results 

Fig. 1(a) presents the high resolution X-ray diffraction (HRXRD) 
pattern of one typical DMO/NSTO thin film, in which pure c-axis 
orientation can be seen for the film. The out-of-plane (OOP) lattice 
parameter of the film was found to be —7.393 A, which is larger than 
that of the bulk counterpart, indicating an in-plane compressive 
strain. Fig. 1(b) shows the rocking curve around the (002) reflection 
for the film, and the full width half maximum (FWHM) is -0.72°. 
The NSTO substrate has a cubic crystalline structure with lattice 
parameter a = 3.905 A, while the orthorhombic DMO has lattice 
parameters a = 5.278 A, b = 5.834 A, and c = 7.378 A 34 . Therefore, 
large lattice misfit is expected in the epitaxial growth. To get detailed 
information of the epitaxial growth, we perform asymmetric recip- 
rocal space mapping (RSM) of the (103) c and (113) c reflections 
where the subscript 'C refers to the pseudo-cubic lattice, and the 
data are shown in Fig. 1(c) and (d), respectively. Identical in-plane 
(IP) lattice parameters of the film with respect to the substrate can be 
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Figure 1 | (a) Symmetric X-ray 6-26 scan of the DMO/NSTO thin film. The inset shows a schematic drawing of the epitaxial growth of DMO on (001) 
NSTO. (b) Rocking curve around the (002) reflection for the film. The FWHM is —0.72°. (c) and (d) present the data of reciprocal space mapping around 
the ( 103) c and ( 1 13) c reflections of the substrates, respectively, (e) 2 X 2 (am 2 atomic force microscopy image of the film, in which islands dominate the 
surface morphology. The R rms of film is estimated to be ~ 1 .6 nm. The inset shows the topography with a smaller scanning size of 0.5 X 0.5 (am 2 , in which 
the island size can be derived to range from 40 nm to 80 nm. 
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seen from the (103) c reflection, indicating the coherent growth of the 
film on the substrate. However, two diffraction spots corresponding 
to (206) o and (026) 0 can be observed in Fig. 1(d), which suggests an 
orthorhombic twin crystalline structure of the film and thus consist- 
ing with the relative large FWHM of (002) reflection. The subscript 
'O' refers to the orthorhombic lattice. Based on the RSM data, a 
sketch of the stacking model of DMO unit cells on NSTO substrate 
is shown in the inset of Fig. 1(a). A similar growth model was prev- 
iously proposed for TMO thin films 35 . Moreover, the atomic force 
microscopy image (Fig. 1 (e)) shows an island surface for the film. 
The root-mean-square roughness (£ rms ) of the film is estimated to be 
1.6 nm, and the island size ranges from 40 nmto80 nm. In addition, 
lattice parameters (and the corresponding variation ratio) of the 
orthorhombic DMO thin films can be estimated from the RSM data 
as, a = 5.288 A (0.2%), b = 5.770 A (- 1.1%), and c = 7.395 A (0.2%). 
Considering the relatively large uniaxial compressive strain along the 
b axis, we expect certain anisotropy in the physical properties of the 
film. 

The T- dependences of magnetization (M) in the ab plane and 
along the c axis are measured under both field cooling (FC) and zero 
field cooling (ZFC) conditions, and the data are shown in Fig. 2(a). 
For both cases, the magnetization evolves smoothly with T until T R — 
9 K, at which distinct anomalies in the M(T) curves are observed and 
can be associated with the long range ordering of Dy moment. 
Usually, the large paramagnetic contribution of Dy moment will 
dominate the magnetic signal, and make the detection of AFM (or 
SSO) phase transition difficult, as reported for bulk DMO single 
crystal 36 . Although identical T R is evidenced from the data in the 
ab plane and along the c axis, the clear difference in shapes of the 
M(T) curves implies the possible magnetic anisotropy within the 
films. In Fig. 2 (b) and (c), the M(H) curves measured under various 
temperatures are plotted. We calculate the demagnetization factor to 
be — 10" 4 in the ab plane and — 1 along the c axis for our sample (area: 
3 mm X 3 mm, thickness: 300 nm), which gives rise to the demag- 
netization field of ~ 1 Oe in the ab plane and — 560 Oe along the c axis 
at H = 7 T 37 . Therefore, the demagnetization field only has very weak 
contribution to our magnetic measurements. Evidently, the magnet- 
ization of ab plane is much larger than that of c axis, which confirms 
the magnetic anisotropy in the films. A close examination of the 
M(H) curves reveals that the maximum magnetization (H = 7 T 
and T < T R ) are M m ab - 10.3 \i B /f.u. for ab plane and M m c - 
1.08 [i B /f.u. for c axis, which gives the magnetic anisotropy as large 
as M m ab /M m c —10. For DMO bulk single crystal, the magnetization 
at H = 7 T and T < T R was reported to be —7.5 [i B /f.u. for b aixs, 
—4.1 [i B /f.u. for a aixs 31 , and —2.5 [i B /f.u. for c axis which is obtained 
by an extrapolation considering the linear M(H) curve 38 , giving rise 
to the magnetic anisotropy to be M m 7M m c — 1.7 and M m VM m c — 3. 
In Ref. 36, the M m b of bulk DMO was reported to be —8.6 \i B /f.u., 
leading to a slightly larger magnetic anisotropy M m VM m c — 3.4. 
Clearly, the magnetic anisotropy of bulk DMO is smaller than that 
of the thin films. Interestingly, a clear metamagnetic transition 
occurs around H = 1.5 T below T R , marked by an arrow in 
Fig. 2(b), which was typically observed in bulk SSO multiferroic 
manganites and associated with the rare-earth moments 8 . It is 
worthy of noting that such metamagnetic phase transition is always 
accompanied by large variation in P due to the one-to-one corres- 
pondence between ferroelectric and magnetic orders 8 . Therefore, we 
explore the possible ferroelectricity and ME coupling in our DMO 
thin films by carrying out detailed dielectric and ferroelectric mea- 
surements as a function of temperature and magnetic field. 

From the T-dependent data of dielectric susceptibility 8 C and 
polarization P c shown in Fig. 3, multiple phase transitions can be 
identified. In Fig. 3(a), a small change in s c arises at T N = 40 K due to 
the ordering of Mn moments. Importantly, a pronounced peak at T c 
= 21 K can be observed in s c (T), and a corresponding emergence of 
P c is identified in the P C (T) plot in Fig. 3(b). According to the phase 
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Figure 2 | (a) Temperature dependence of magnetization measured in the 
ab plane and along the c axis. Measured independence of magnetization in 
the ab plane (b) and along the c axis (c) at various temperatures. The arrow 
in (b) marks the metamagnetic transition. 

diagram of orthorhombic perovskite manganites, two magnetic con- 
figurations, i.e. the SSO and E-type AFM orders, can break the spatial 
inversion symmetry and thus produce a spontaneous polarization 22 , 
while the latter has a spontaneous polarization along the a -axis 39 ' 40 . 
Moreover, a recent first-principles calculation predicted that the E- 
AFM multiferroic phase can only be induced by large tensile strain 
(5%) 41 . Therefore, the polarization P c of the compressively strained 
thin films should origin from the SSO via the DM interaction, similar 
to the bulk counterpart 9 " 11 . The SSO induced ferroelectric state is 
further attested by the well-shaped polarization -electric field (P-E) 
hysteresis loop, as shown in Fig. 3(c). Generally, in SSO multiferroics, 
the reversal of polarization from +P C to — P c corresponds to the 
reversal of spin chirality (Q = 5 2 - X Sj) from clockwise Q_ a (be-) 
to counterclockwise Q a (bc+) (or vice versa), marking the close link 
between magnetic and ferroelectric orders 42 . 

Further cooling the film down to T — 12 K, both s c (T) and P C (T) 
(red circles in Fig. 3(b)) abruptly decrease, indicating a possible trans- 
ition between different magnetic configurations. Actually, similar 
features were also observed in bulk DMO, and a SSO to commen- 
surate (CM) AFM phase transition of Dy moment was considered as 
the origin 30 " 33 . In details, Dy moments form a SSO phase with the 
same propagation vector as the Mn magnetic lattice (t Dy = T Mn = 
0.385b) at 15 K (<T C - 18 K), due to the Dy-Mn spin interaction, 
and then make an additional contribution to P c 30 . Further cooling 
causes a suppression of the Mn-promoted SSO of the Dy moment, 
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Figure 3 | Temperature dependence of e c (a) and P c (b). The arrows in (a) indicate the directions of the heating or cooling processes. Red circles and black 
squares in (b) correspond to the results of different poling processes. In (c), ferroelectric hysteresis loop measured at T = 9 K without magnetic field 
during the cooling process. In the upper of (c) , a sketch of the clockwise ( Q_ J and counterclockwise ( QJ spiral spin orders is presented, corresponding to 
the +P c and — P a respectively. 



and a CM/AFM configuration with T Dy = 0.5b arises because of the 
dominance of the Dy-Dy interaction over the Dy-Mn spin coupling 
at T < T R 31 " 33 . As a consequence, abrupt decrease in P c and s c occurs, 
as observed in the P C (T) and s c (T) data of the film (Fig. 3(a) and (b)). 
Moreover, microscopic evidences revealed significant thermal hyster- 
esis for the first order SSO -CM/AFM phase transition of Dy 33 , which 
suggests a strong thermal history dependence of the microscopic 
process related to the emergence of FE polarization. To confirm this 
point, we poled the films using the same electric field down to 8 K 
(close to T R ), and then the polarization was collected by integrating 
the pyroelectric current. Remarkably different P c data evolving with T 
for the two measuring processes were indeed obtained as shown in 
Fig. 3(b), indicating the important role of Dy moment to the 
improper ferroelectricity in the films. 

To further explore the coupling between magnetic and ferroelec- 
tric orders, which in turn would be helpful for us to comprehensively 
understand the physics of multiferroicity in DMO thin films, we 



measure the T- dependent e c and P c under various magnetic fields. 
As displayed in Fig. 4, rather different magnetic responses of s c (T) 
and P C (T) can be identified when H is applied along the c-axis or in 
the ab plane. For Hllc axis, no substantial H dependence in P c is seen 
although T c is slightly suppressed with increasing H. For HI lab 
plane, remarkable enhancement in P c is observed, and a maximum 
value of P C (2T)AP C (0T) — 9 is estimated at T = 2 K. Correspondingly, 
the abrupt decrease in s c below 12 K, induced by the long range 
ordering of Dy moment, is significantly suppressed by HI lab plane, 
indicating the essential role of Dy moment in mediating the FE 
polarization in the films. At H > 5 T, the FE phase transition in 
s c (T) becomes dispersive, and as a consequence P c is largely sup- 
pressed even in the temperature region T > T R . It implies possible 
P- flopping due to the rotation of the frc-spiral to the afr-spiral caused 
by HI lab plane. 

In Fig. 5, the H dependence of pyroelectric current J, P c and s c (H)l 
e c (0) at T = 2 K are displayed. Distinct current peaks and valleys in 
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1(H) signal the change (for example flopping or onset) in P c related to 
the variation of spin configuration induced by the magnetic field. 
Accordingly, we can roughly partition the data into three regions, 
namely, phase I (H < 1.5 T), phase II (1.5 T < H < 4.5 T), and phase 
III (H > 4.5 T). In phase I, the CM-AFM phase of Dy moment 
coexists with the SSO phase of Mn, and P c totally originates from 
the spirally ordered Mn sublattice, following the DM mechanism. 
Around H ~ 1.5 T, coinciding with the metamagnetic phase trans- 
ition of Dy in M(H) (Fig. 2(b)), the integrated FE polarization P c (in 
Fig. 5(b)) exhibits an abrupt increase, which is accompanied by 
anomaly in s c (H) (in Fig. 5(c)). It suggests the direct effect of the H 
induced new magnetic orders of Dy on the ferroelectricity, while 
microscopic evidence is needed to determine whether it is the spiral 
spin order similar to the case in bulk DMO. In phase II, multicom- 
ponent magnetic structures coexist with each other, and P c reaches a 
plateau. As H > 4.5 T in phase III, P c is gradually suppressed and 
hysteresis anomaly can be observed in £ C (H), likely indicating the 
magnetic field induced rotation of the frc-spiral (P c ) to the afr-spiral 

(Pal 

A key feature in Fig. 5 is the emergence of the new spin ordering 
phase of Dy around 1.5 T, which is accompanied by distinct enhance- 
ment in P c . This feature implies dramatic change in the multiferroic 
domain structure, and thus leading us to expect a different ferroelec- 
tric switching behavior along with the phase transition of Dy. 
Fig. 6(a) shows the P-E loops measured under T = 5 K and various 
magnetic fields, in which clear difference in the FE coercive field (E c ) 
can be observed besides the different polarizations upon increasing 
HI lab plane. In Fig. 6(b), the P c and E c is plotted as a function of H. 
As expected, the P c shows a sudden increase around H = 1.5 T 
(corresponds to the metamagnetic phase transition of Dy moment), 
which is well consistent with the pyroelectric current data shown in 
Fig. 5. Another striking feature is the dramatic increase in E c around 
H = 1.5 T, which implies a higher energy cost in the ferroelectric 
switching (or chirality Q switching). A sketch of energy diagram on 
the FE switching process is plotted in Fig. 6(c). It is clear that the 
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Figure 5 | Isothermal pyroelectric current I c (a), integrated P c (b), and s c 
(c) as a function of magnetic field HI lab plane measured at T = 2 K. The 

orange dots are the date obtained from the P c ( T) measurements. 
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Figure 6 | (a) Ferroelectric hysteresis loops measured under various 
magnetic field at T = 5 K. (b) independence of the switching polarization 
P c (open squares) and ferroelectric coercive field E c (open circles) 
measured using the PUND method. The orange hatched region denotes 
the magnetic field driven phase transition of Dy moment, (c) A sketch of 
the energy diagram on the ferroelectric switching process. 

reversal of chirality Q from Q a (Q- a ) to Q- a (Qa) should overcome 
the barrier coming from the energy difference AE between Q a /- a (be 
spiral) and Q c /- C (ab spiral) 43 . Therefore, the larger E c suggests a 
higher energy barrier between the ab and be spirals (indicated by red 
dot line in Fig. 6(c)) caused by the emergence of Dy spin order. 

Discussion 

To this stage, we have demonstrated the magnetic ferroelectricity and 
ME coupling in DMO thin films. In qualitative sense, the DMO films 
do show similar physical properties to the bulk counterpart, includ- 
ing the anisotropic magnetism, the SSO phase induced ferroelectri- 
city, and the magnetic field induced polarization enhancement 8 ' 22 . 
These similarities imply the similar physics of multiferroicity for 
both thin film and bulk DMO. The present work marks a successful 
copy of the bulk properties to the thin films for SSO multiferroic rare- 
earth manganites. In fact, contrast to the case of DMO thin films 
shown here, so far no ferroelectric polarization has been reported in 
the TMO thin films 13 " 18 . A possible reason for this difference is that 
the DMO locates deeply inside the SSO phase, while the TMO locates 
in the vicinity of the A-type AFM phase 22 . Therefore, the SSO phase 
in TMO thin films can be sensitively influenced by perturbations 



SCIENTIFIC REPORTS | 3 : 3374 | DOI: 1 0.1 038/srep03374 



5 



such as strain, while the SSO phase in DMO thin films shows rela- 
tively high stability against the strain. 

Even though the substantial similarity in those multiferroic beha- 
viors between the DMO bulk and thin film, one still can identify 
several quantitative differences in these behaviors between them. 
First, for DMO thin films, the in-plane and out-of-plane magnetiza- 
tions are different from those of the bulk as aforementioned, giving 
rise to a larger magnetic anisotropy M m ab IM m c ~ 10 in the thin films. 
Based on the phase diagram and crystalline structure of perovskite 
manganites, we note that the alb ratio exhibits the same trend of 
evolution with the GdFe0 3 -type distortion and the Jahn-Teller (JT) 
effect as the A- site rare -earth element is changed from La to Er, which 
means the larger alb value the smaller GdFe0 3 -type distortion and JT 
effect 34 ' 44 . Therefore, the larger alb ~ 0.916 in the DMO thin films 
implies the partial suppression of the JT and GdFe0 3 -type distortion 
by the IP compressive strain, as compared with the bulk DMO in 
which the alb ratio is 0.905 34 . Consistently, the larger alb ratio 
accompanied by the suppression of the JT distortion was also reported 
for bulk DMO by synchrotron XRD experiments carried out under 
static pressure 45 , which is analogous to the compressive strain in our 
thin films. Moreover, it is well known that the large JT and GdFe0 3 - 
type distortion are responsible for the SSO phase in multiferroic 
manganites 9 " 12 . Thus, the compressive strain suppressed JT and 
GdFe0 3 -type distortion will weaken the spin-spirals (also the DM 
interaction) and give rise to a reduced P c (—31 |iC/m 2 at 2 K and 
H = 0) in the thin films. From another perspective, according to the 
model of Mochizuki et a/. 10 " 12 , which can theoretically reproduce the 
phase diagram of orthorhombic multiferroic manganites, the DM 
interaction coexists and competes with the single-ion anisotropy in 
a way that the frc-spiral is favored by the former while the afr-spiral is 
stabilized by the latter. Because of the weakening of DM interaction by 
strain in the films, more pronounced magnetic anisotropy is thus 
anticipated, in comparison with the bulk counterpart. Indeed, the 
magnetism of both bulk and thin film manganites was revealed to 
well scale with the alb value 46 , which further confirms our explanation 
on the strain mediated magnetic anisotropy in the DMO thin films. 

Besides the magnetic anisotropy, an enhancement in P c (P C (2T)/ 
P c (0T) ~ 9 at T = 2 K) is observed in the DMO thin films, as shown 
in Fig. 4(d), which is larger than that (P C (2T)/P C (0T) ~ 3 at T = 2 K) 
in the bulk counterpart 22 . On one hand, because of the large magnetic 
anisotropy in the DMO thin films, the possible Dy-spiral phase 
indeed deserves further microscopic characterization. On the other 
hand, the above outlined experimental results reveal direct effect of 
Dy moment on the ferroelectric polarization, which signals the Dy- 
Mn spin interaction and thus leads us to argue a coherence of spin 
configurations for Dy and Mn. In fact, the coherence of R and Mn 
was extensively identified in multiferroic RMn0 3 (R = Tb, Dy, Ho) 
due to the R-Mn spin interaction 31,33 ' 47,48 . Such magnetic struc- 
ture meets well the configuration of symmetric magnetostriction 
effect 49 " 51 , that is, each Dy spin chain is sandwiched by two antifer- 
romagnetically coupled Mn spiral spin chains with the same pro- 
pagation vector. The stacking of the magnetic orders is along the c 
axis, and the Dy orders would displace cooperatively toward (away) 
the Mn orders with antiparallel (parallel) spins via the symmetric 
magnetostriction effect and, as a consequence, polarization along the 
c axis would be generated, as sketched in Fig. 7(a). Actually, the 
symmetric magnetostriction effect was proposed as the origin of 
the polarization enhancement in bulk DMO before the Dy-spiral 
phase was demonstrated experimentally 52 . Moreover, the symmetric 
magnetostriction effect between coherent Dy and Mn magnetic 
orders would also be expected in bulk DMO which has Dy-spiral 
phase exhibiting the same propagation vector with Mn SSO, as 
shown in Fig. 7(b). However, the different magnetic anisotropics 
for thin film and bulk DMO would induce different magnetostriction 
effects on the polarization, which is the possible origin of the H- 
induced enhancement in P c shown in the present work. 




M C 



O Dy o Mn O 

> 

Figure 7 | Sketch of the coherence of magnetic configurations for Dy and 
Mn spins, (a) the magnetic order of Dy just has the same propagation 
vector as Mn-spiral phase, (b) Dy has a spiral spin ordering as Mn spins. 
The both configurations meet the magnetic structure of symmetric 
magnetostriction effect, and thus possibly leading polarization (P s ) along c 
axis due to the uniform displacement of Dy. 

In summary, multiferroic DyMn0 3 thin films are epitaxially 
grown on Nb-SrTi0 3 (001) substrates, and intensively characterized. 
Our results reveal intriguing multiferroicity in DMO thin film, 
including i) large magnetic anisotropy M m ab IM m c —10 below T R , 
which is induced by the compressive strain suppressed GdFe0 3 -type 
distortion and Jahn-Teller effect in the thin films; ii) large enhance- 
ment in P c (—800% at T = 2 K) and FE coercive field E c , which was 
proposed to be understood by interacting magnetic structures. Our 
work shows a successful copy of the spiral spin ordering induced 
multiferroicity from the bulk to thin films, which opens the possibil- 
ity to integrate such materials into functional devices. 

Methods 

Pure phase DMO thin films with thickness t = 300 nm were epitaxially deposited on 
(001) Nd-SrTi0 3 (NSTO) substrates using pulsed laser deposition (PLD) at 800°C in 
0.2 mbar oxygen. KrF excimer laser (X = 248 nm) was used with an energy density of 
~ 1 J/cm 2 . After the growth, the samples were cooled to room temperature in 400 
mbar oxygen ambient at a rate of 5°C/minute. Conducting NSTO substrates were 
used as both growth templates and bottom electrodes. Detailed crystalline structure 
investigations, including the normal 6-29 scan, rocking curve, and the asymmetric 
reciprocal space mapping, were carried out using a high resolution X-ray diffraction 
system (Panalytical X-pert Pro). The surface morphology of the film was probed 
using tapping mode atomic force microscope (Digital Instrument). 

Magnetization as a function of T and magnetic field (H) was measured using a 
Magnetic Properties Measurement System (MPMS, Quantum Design). The cooling 
and measuring field during the M(T) measurements was fixed at 0.1 T. To obtain the 
FE polarization along the c axis P c of the films, square Au top electrodes of 300 um in 
diameter were prepared on the film surface, and the pyroelectric current method was 
used to probe the polarization. For each measurement, the sample was first cooled 
from 60 to 2 K under a poling electric field of 100 kV/cm. Then, the poling electric 
field was turned off and the sample was short-circuited for a sufficient period of time 
to release the charges accumulated on the sample. The pyroelectric current was 
measured using a Keithley 6514A electrometer during the heating process with a 
constant ramping rate (2 K/min to 5 K/min). For the pyroelectric current measure- 
ments under a magnetic field, H was applied during the poling and the subsequent 
measuring processes. The ferroelectric hysteresis loops were measured using the 
positive-up-negative-down (PUND) method (TF Analyzer 2000 from AixACCT Co), 
and the frequency was set as 100 Hz. The dielectric susceptibility s c of the sample with 
the same device structure was measured using a HP4294A impedance analyzer. The 
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cryogenic environment and the magnetic field during the pyroelectric current and 
dielectric susceptibility measurements were provided by a Physical Properties 
Measurement System (PPMS, Quantum Design). 
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